A NUMERICAL SCHEME FOR IMPACT PROBLEMS 

LAETITIA PAOLI * AND MICHELLE SCHATZMAN t 

Abstract. We consider a mechanical system with impact and n degrees of freedom, written 

in generalized coordinates. The system is not necessarily Lagrangian. The representative point is 

^^ ^ subject to a constraint: it must stay inside a closed set K with boundary of class C^. We assume 

0\ 1 that, at impact, the tangential component of the impulsion is conserved, while its normal coordinate 

Qs ■ is reflected and multiplied by a given coefficient of restitution e £ [0, 1]: the mechanically relevant 

notion of orthogonality is defined in terms of the local metric for the impulsions (local cotangent 

metric). We define a numerical scheme which enables us to approximate the solutions of the Cauchy 

problem: this is an ad hoc scheme which does not require a systematic search for the times of impact. 

O, We prove the convergence of this numerical scheme to a solution, which yields also an existence 

■ result. Without any a priori estimates, the convergence and the existence are local; with some a 

priori estimates, the convergence and the existence are proved on intervals depending exclusively on 
these estimates. The technique of proof uses a localization of the scheme close to the boundary of 
K; this idea is classical for a differential system studied in the framework of ffows of a vector ffeld; it 
^^ , is much more difficult to implement here, because ffnite differences schemes are only approximately 

^^^ . local: straightening the boundary creates quadratic terms which cause all the difficulties of the proof. 
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1. Introduction. We study in this article a numerical approximation of dynam- 
^ ' ics with impact with a finite number of degrees of freedom and a smooth constraint. 

CO . The set of constraints is denoted K and satisfies the following assumptions: 

^^ , K is a closed subset of R'^ with non empty interior; (1.1a) 

{the boundary dK of K is an embedded sub-manifold , , , 

, ^ (1.1b 

of class C3 of M''; ^ ' 

^ i K lies on only one side of dK. (1-lc) 

It is possible to find a function </> of class C^ such that 

K = {ueR'^ : 0(w) > 0} 

'V^ . and the differential d(f> does not vanish on dK = {w G M'' : 0(w) = O}. 

Cd ' Let / be a continuous function from [0, T] x M*^ x R'^ to M'' which is locally Lipschitz 

continuous with respect to its last two arguments, and let M{u) be the mass matrix: 
u H^ M{u) is a mapping of class C^ from W^ to the set of symmetric positive definite 
matrices. 

The free dynamics of the system are written in generalized coordinates as 

M{u)u^ f{-,u,p), p^M{u)u. (1.2) 
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This system is more general than the system obtained in Lagrangian mechanics, 
since we want to include possible dissipative terms in the dynamics of the problem 
under discussion. 

Let us give the few geometric notations which are absolutely necessary here, since 
we use a Riemannian metric; the cotangent bundle T*M.'^ is identified to M.'^ x M.'^, and 
its elements are denoted as pairs (u,^); at each point u of M."^ the metric tensor for 
tangent vectors is defined by the matrix M{u), and the metric tensor for cotangent 
vectors is defined by the matrix M{u)~^. The scalar product of two vectors x and y 
in the tangent space at u is denoted by (a;,?/}„; coordinate-wise it can be expressed 
as x^ M{u)y where x and y are column vectors. The scalar product of two vectors 
^ and rj in the cotangent space at u is denoted by (Cj^}u ^'^'^ coordinate-wise it is 
equal to ^'^ M{u)~^rj. The corresponding norms of vectors and covectors are denoted 
respectively by |a;|tj and |^|* . 

Therefore, a cotangent vector (u, f ) belonging to T*W^ is orthogonal to the cotan- 
gent vector [u.-q) iff (C,??)^ vanishes. 

With these notations, if the velocity of the system is u, the generalized impulsion 
is M{u)u = p and (m,p) belongs to the cotangent space T*W^. Whenever we take the 
orthogonal of a vector or a vector subspace of the tangent or the cotangent space at 
u, we always use the relevant metric tensor; therefore it is important to know which 
of the vectors under consideration are cotangent and which are tangent. Of course, 
all the differential forms are cotangent vectors. 

Let us describe now the system satisfied by the problem with impact: we re- 
place dT]) by 

M(u)u = ^i + f{-,u,p), (1.3) 

and since we cannot expect to have global solutions in general, /.t is an unknown 
measure on [fo, ^o+t^] with values in W^ which describes the reaction of the constraints: 
// has the following properties: if rf^ denotes the differential of 0, then 

supp(/x) C {t e [to, to + r] : 0(w(t)) = 0}, (1.4a) 

fi^\d(l){u), (1.4b) 

A > almost everywhere on [to, to + t]. (l-4c) 



We require the following functional properties for u: 

{M is a continuous function taking its values in K 
for aU t G [to,to + T], 

u is of bounded variation over [to, to +t]. (1.5b) 



(1.5a) 



If it is of bounded variation, p is also of bounded variation. Assume that u{t) 
belongs to 97^; we decompose p(t — 0) andp(t-l-O) on K(i0(w(t))©(i(/)(ii(t))^; here the 
_L sign means the orthogonality with respect to the local cotangent metric. We inte- 



grate (1.2) on a small neighborhood of t, relation ( 1.4b ) implies that the component 
oi p{t — 0) on d(l>{u{t))-^ is conserved. 

Therefore, we have to make a supplementary assumption in order to have a com- 
plete description of the impact; we choose a constitutive law of the impact using a 
coefficient of restitution: thus we will assume that there exists e G [0, 1] such that the 
component of p(t + 0) along M.d(t){u) is equal to — e times the component of p{t — 0) 
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on Rd(p{u). In other words, we have 

Pit + 0) = Pit - 0) - (1 + e) fU ^uit)). (1.6) 

W(w(i)),#(w(0))^(t) 

The set of admissible initial data D will be 

P={(to,uo,Po) e [0,r)xif xR'^: 
if Wo e dK, then (po, ^^('"o))*^ > o|. 



(1.7) 



This choice is equivalent to the convention that there is no impact at the initial time 

to. 

Given initial conditions (io,wo,po) G 15, we require that the following Cauchy 
data be satisfied: 



and 



uita) = uq, 



p{to) =pa. 



(1.8) 



(1.9) 



For all initial data (to, ua,Po) G © we will obtain the existence of a local solution 



to (1.2), (1.4a), (1.41), (1.4c) and (l.E) belonging to the functional class defined 
by (1.5a) and (1.5b) and satisfying the initial conditions (l.S) and (1.9). 

The existence of this local solution is obtained by defining a numerical scheme, 
whose convergence will be shown in appropriate functional spaces; the limit of the 
approximation will be a solution of our problem. 

The distance on M'' is defined with the help of the Riemannian metric: if s ^-> u{s) 
is a C^ mapping from [a, b] to M"*, the Riemannian length of the image of u is 



i{u) 



l"(s)L(s) ds. 



This curve length is invariant by a diffeomorphic change of parameter. Therefore, we 
may assume that a = and b — 1. The distance from x to y is the lower bound of 
the length of the curves from x to y, or in other words: 



dist(x,2/) =inf{£(u) : ueC\[0,l]), u(0) = 



41) = y}- 



It is classical that the lower bound is attained on the geodesies for the given Rieman- 
nian metric; it is also known that for each point x there exists r > such that if 
dist(a;, y) < r there is only one geodesic from x to y. 

We denote by dist(a;, E) the Riemannian distance of a point a; to a set E. 

Under assumptions ( |1 . l| ) , a projection on dK can be defined uniquely on an 
appropriate neighborhood of dK; more precisely, for all compact C C dK, there 
exists a neighborhood of C on which the projection Pqk is uniquely defined, and 
there exists a unique geodesic joining a point of this neighborhood to its projection. 
This projection Pqx is characterized by the relation 



Vy G dK, distiPsKX, x) < dist(2/, x). 



(1.10) 
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This projection is of class C^. 

For all X in dK, denote by N{x) the interior unit normal vector: this means that 
\N{x)\x is equal to 1 and that it is orthogonal to the tangent space at Pokx with 
respect to the scalar product in the tangent space, i.e. for all y such that d(f){x)y 
vanishes, {y,N{x))x = 0. The smoothness of dK implies that the mapping z i— > N{z) 
is of class C^. 

When the geodesic from x to Pokx is unique it is tangent at Pokx to N{PdKx). 

Starting from this projection on dK, we can define a projection on K as follows: 
for each compact C included in K, there exists a relatively compact neighborhood U 
of C on which Pk is defined by 

Pk{x) = < . (1.11) 

I X otherwise. 

The reader will check that Pk is Lipschitz continuous over U and that Pkx realizes 
the minimum of the distance from x to K. 

Given two positive numbers h* <1 and T, assume that i^ is a continuous function 
from [0,r] X R'^ X K'* X K'^ X [Q,h*] to W^, which is locally Lipschitz continuous 
with respect to its second, third and fourth arguments; assume moreover that F is 
consistent with /, i.e. that for all t £ [0,T], for all u and v in M'* 

F{t, u, u, V, 0) = M{uy^f{t, u, M{u)v). (1.12) 



We approximate the solution of (|l|), ( pTIa] ), (|L4|), ( pTI^) , ( pTSa] ), ( pTSb] ), (^), 



(1.9) by the following numerical scheme: the initial values U'^ and U^ are given by 



the initial position 

C/° = Mo, (1.13) 

and the position at the first time step 

U^ ^uo + hM{uor^Po + hz{h), (1.14) 

where z{h) tends to as ft. tends to 0. 

We will use systematically henceforth the notation 

tm = to + mh. (1-15) 

Given U"^~^ and [/'", [/™+i is defined by the relations 

t/"- = ^eU-^ + (1 + e)P. (^^U--i^-e)U--^ + ^^'F-^^ (l.,6) 

and 

/ rrm+l _ Tjm-l \ 

F"' = F(t^,U"',U"'-\ ,h\ (1.17) 

provided that C/™"'"^ is unique in a neighborhood of t/™. 

A commentary on the construction of this scheme from the point of view of 
convex analysis will be useful here. We refer to the book of Rockafellar [g^] for more 
information on the basic ideas in convex analysis to be used below. 
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Let us assume provisionally that the set of constraints K is convex and that the 
mass matrix is equal to the identity matrix on W^ . Then the Riemannian structure 
of W^ is simply its Euclidean structure. 

Recall that the indicator function iI)k of a closed convex set K is defined by 



iPk{x) = 







if a; e K, 
-oo otherwise, 



(1.18) 



and its sub-differential dipK is a function from K to the set of closed convex sets given 

by 



dipK{x) 



{0} 



-N{x) 



if X e \nt{K), 
if X e dK. 



For all A > 0, the multivalued equation 

X + XdtpK^x) 3 f 
has a unique solution given by 

X = PKif), 



(1.19) 



(1.20) 



(1.21) 



where Pk is the usual projection on the closed convex set K in Euclidean M''. 

In the annoucement ||2^, where we assumed that the set of constraints K was 
convex and the geometry was Euclidean, we had defined the numerical scheme by the 
multivalued equation 



[/™+i - 2C/" 



[/" 



/l2 



dtPK I :-^ I 3 F" 



(1.22) 



We may rewrite ( 1.22 ) as 



-dipK 



Jjm+1 j^ e^m-1 



1 +e 1 +e 

2C/™- (l-e)C/™~i +h^F" 



1 



(1.23) 



1 



e 



which reduces, thanks to (1.20) and (1.21) to relation (1.16). 

If we generalize (1.22) to a non convex K with a general mass matrix, we cannot 
use the apparatus of convex analysis, and there is no good reason to use the even 
more technical apparatus of non-convex analysis a la Clarke: this theory is useful 
when the corners of K are not convex; in the mechanical setting, corners are convex, 
since they appear as the intersection of smooth sets of constraints. Here, the problem 
is even simpler because we do not have any corners. 

The boundary dK is smooth, and as we expect that for small h, the [/'"'s will 
stay close to K, we still have a projection of (2t/™ - (1 - e)U"'-'^ + /i2f"')/(1 + e) 
on K, and thus we start from (1.16) to define the numerical scheme. 

The original definition reappears as follows: define 



W 



2C/™ - (1 - e)[/™-i + h^F" 
TT~e 



(1.24) 



L. PAOLI AND M. SCHATZMAN 



that will be used in many places in the upcoming proofs. With this definition, ( f .fq ) 
is rewritten as 



U 



rn+l 



-et/'"-i + (l + e)P;f(W^"). 



Hence, if we define 



we find that 



U 



711 -\-l 



eU 



m—1 



l + e 



Pk{W"'). 



If we subtract ( 1.25 ) from ( 1.24 ) , we can see that 

jjrn+l _ 2UTn _^_ jjni-1 ^ 



/l2 



h^ 



(VK" - Z™) = F" 



(1.25) 



(1.26) 



(1.27) 



which reduces to (1.22) in the convex case with a trivial mass matrix. 

Another way of writing (1.27) is to define the discrete velocity y™ by 



y^ 



U 



m+l 



jjn 



Then, ( |l.27 ) can be rewritten as 



yni _ yr. 



hF""^ 



(l + e)(Z'"-Ty") 



(1.28) 



(1.29) 



A strict contraction argument in M"^ gives the existence of a unique C/™ for small 
values of m and h. As the projection on K is uniquely defined only in a neighborhood 
of K, and is only Lipschitz continuous, the iteration of a fixed point argument might 
request smaller and smaller bounds on the time step h, and there is no guarantee that 
we could integrate numerically on a time interval bounded from below, for any initial 
time step size. 

It should be noted that this difficulty is specific to the non convex case. 

Let us outline now the structure of the article and of the proofs. In the one- 
dimensional case, the main estimates are given by lemma 2.1 , in section 0. In sectionH, 
we will straighten the boundary, a natural geometrical idea. 

While the syst em (p^)-(|1.6D is nicely transformed under a diffeomorphism, the 
numerical scheme ( LlBj ), (1.14), (1.16) and (1.17) does not behave well under diffeo- 
morphism. The reason is that a numerical scheme is not a local object: when we define 
a discrete velocity by subtracting C/™ from U™~^^, we use locally a vector structure 
which is not intrinsic from the point of view of differential geometry. In particular, if 
we apply a diffeomorphism to the numerical scheme, we will find another numerical 
scheme which will look much more complicated than the previous one, since it will 
contain a number of small term which show the lack of an intrinsic description of the 
scheme. After a very technical proof, we find two constants C3 and r such that for 
initial data in a compact subset of the admissible set, and for all small enough h and 
all m < r/h, the discrete velocity is bounded: 



sup IF™! <C3. 
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Since uniqueness is not true in general , |£9| , and hypotheses of analyticity are 
often but not always used for the proof of uniqueness [2^ , [g^ i Q ' Hi ' ' the proof 
of convergence of the numerical approximation is delicate also for this reason. 

However, there is a bonus: all the effort made to prove the local convergence of 
the numerical scheme provides us with a local existence proof for our problem. In 
sections §, ^ ^ and 0, we prove estimates on the discrete acceleration, we establish 
the variational properties of the limit of the numerical scheme, and we study the 
transmission of energy at impact, as well as the passage to the limit for the initial 
conditions. All these results are obtained under the assumption that on a certain time 
interval starting at to: the discrete velocity is bounded independently of the time step. 

As a pre li mina ry to the global existence proof, we give a priori estimates on 



problem (1.3)-(1.£) in section R which is completely independent from the remainder 
of the article. 

In section 0, we establish a very weak semi-continuity for the supremum of the 
local norm of the discrete velocities; this result enables us to obtain a global existence 
and convergence theorem. 

This article is of a theoretical nature: the existence result obtained here is a 
generalization of [|9), |], |§, (l|, |§. 

The numerical scheme analyzed here has been implemented in the case of a trivial 
mass matrix in [n9|, 123], 118], |23]. In all these articles, we compared the performances 
of this scheme with those of a method based on the detection of impact. When the 
impact times are isolated, the scheme by detection of impacts is more precise than 
the present scheme. As soon as the restitution coefficient is strictly less than one, 
we find systematically non-isolated impact times. In all cases, the present scheme is 
substantially faster. Since the phenomena that we want to approximate are highly 
nonlinear and often very sensitive to the initial data, the issue of precision is not 
necessarily crucial. Our numerical experiments show that the performance of the 
present numerical scheme is quite satisfactory from the point of view of qualitative 
conclusions. 

The case of a non-trivial mass matrix, and a stiff system, indeed the case of the 
discretization of a beam has been adressed in [E4| . 

Let us remark that many articles have been devoted to the problem treated here, 
under the assumption of anelasticity, i.e. a situation where the normal component 
of the impulsion vanishes after the impact; Moreau applied Gauss' principle of least 
constraint to unilateral problems in order to justify his choice of anelastic impact |l3], 
which eventually led him to sweeping processes U^ , followed by |13| , p4[ ; dry friction 



enters in Moreau's work as ||16|; frictionless anelastic impact starts as 17 , and the 
mathematical theory is tackled by M. Monteiro-Marques in a series of articles: his 
main contributions are n^ for the general theory of differential inclusions, um for 
one-dimensional dynamics with friction, [p| which adds percussion to the previous 
framework; this work is improved as |9[|, where dynamics of n particles on a plane 
with normal friction are considered. The discretization approach has been taken up by 
Monteiro-Marques and Kuntze in M , but most significantly by Stewart and Trinkle: 
they use that approach in [ 31 1 , |3^ and [g^ ; the real coronation is the beautiful and 
difficult article of Stewart |32], which concludes the study of dynamics with friction 
and anelastic impact for a finite number of degrees of freedom, and one constraint, 
and still important results in the multiple constraint case. 

The philosophy of this long list of works is somewhat different from ours: we feel 
that not all impacts are anelastic, and we were originally motivated by continuous 
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media; thus, we wanted to develop methods which work weU for stiff systems of 
ordinary differential equations. From this point of view, any method which has to 
calculate with some precision the impact times is doomed to failure. On the other 
hand, the precision of the method presented here needs improvement, and globally, it 
would make sense to agree on benchmarks which would enable the end-user to decide 
between different numerical methods. 

2. The heart of the estimates. In the one-dimensional case, the main estimate 
on the numerical scheme is described in the following lemma; we recall the definition 

r^ = max(r, 0). 

Lemma 2.1. Let the real-valued sequence (y™), satisfy the following recurrence 
relation for all m > 1 : 

y™+i = -ey"-^ + (2y'" - (1 - e)y'"-i)"*' + h^X"". (2.1) 

Then, for all m>2, the discrete velocity 

<' = (2/"+'-y™)//i (2.2) 

satisfies the estimate 

|ry™| <max(|r/™-i| ,e\r^"'-'^\) + h\\"'\ + h\y"-^\ . (2.3) 

Proof. Assume first that 2y™ — (1 — e)w™~^ is non negative, and substitute 

y-m+l ^ y-,n _^ i^^m.^ym-1 = ym _ /j^m-l ^^^^ (|^). ^g obtain 



so that 

h"! < h^'M +/t|A™|. (2.4) 



Assume now that 2j/™ — (1 — e)?/™ ^ is strictly negative. On one hand, (2.1) implies 
the relation 

h 
the assumption on the sign of 2^™ — (1 — e)j/™~^ is equivalent to 

^i^<-(l-e)<-, 

and therefore 

7?" > 7?"-! + /lA". (2.5) 

On the other hand, we subtract from the relation 



SCHEME FOR IMPACT 



the inequality implied by ( p.l| ) with m substituted by tti — 1: 



y" + ey 



™-2 > h'^X"'-^ 



and we infer that 



< - 



erj 



ni-2 



+ /i(A'"-A'"-i). 



(2.6) 



When we summarize (2.4), (2.5) and (2.6), we find (2.3). D 

Later on, we will give a d-dimensional version of (2.3), where the main difference 
is due to geometric effects: there will be a term resembling A™, and the game will be 
to prove a bound on this term. 

3. Existence of (J7™)o<m<[r//ij for some r > 0. We use systematically the 
floor and ceiling notations: when r is a real number, the floor \r\ of r is the largest 
integer at most equal to r, and the ceiling [r] is the smallest integer at least equal to 
r. 

The main result of this section is the existence of a number r > such that for 



all small enough h and all m < [r//ij there exists indeed a discrete solution of ( 1.16| ) 
and ( 1.17 ), whose discrete velocity is bounded independently of h. In fact, we prove 
a stronger result: provided that the first two discrete velocities are bounded, we find 
a uniform lower bound on r when the initial position belongs to a compact subset of 
K. 

We prove first the existence of C/^ under appropriate assumptions on U^ and U^. 
This proof decomposes in two lemmas: the first lemma is strictly an initial condition 
statement, in which no uniformity with respect to initial conditions can be obtained. 
The second one will be used in the foregoing induction proofs. 



Lemma 3.1. For all {tQ,uo,M{uo)vo) e D, f or all U^ satisfying (|l.l4| ), and for 
all small enough h, there exists a solution U'^ of ( 1.16 ) for m = 2 satisfying 



\U'-U\<2\vo\uoh. 



Proof. Let r > be such that Pk is Lipschitz continuous on 

Buoiuo,r) = {ueW^ : |m-woL, <r}. 

Define Ci by 

Ci^nmx{\F{t,u,u',0,h)l^ : ie[0,T], \u - uo\^^ < r, 
\u'-uo\^^<r,he[0,h*]}, 

and let L be the Lipschitz constant defined by 

~ f\F{t,u,u',v,h)~F{t,u,u',v',h)l^ 

L = sup -— — : \u - uoL„ < r, 

I"' - uo\^^ < r, \v\^^ < 2 |i;oL„ + 1, \v'l^ < 2 |^,oL„ + 1, 
vj^v', he [0,/i*]j. 
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Finally, let 7 be the Lipschitz constant of Pk defined by 
(\Pku-Pku'\ 



7 — sup ■ 



\u — u 



\u - uq\ <r, \u- uqL <r, u^u 



There exists a function z{t) which is bounded in a neighborhood of such that for 
small positive values of t: 



Pk{uo + tvo) — uo+ tvQ + t^z{t)] 



(3.1) 



indeed, if wq vanishes, or if uq belongs to int(-fir), or if uq belongs to dK and the 
scalar product {vq, N{uo))uo is strictly positive, z va nishes; if uq belongs to dK and 
(wo, N{uq))uo vanishes, while vq does not vanish, (3.1) is a consequence of the smooth- 
ness of Pqk in a neighborhood of uq : for the values of t for which uq + tvo belongs to 
K, z vanishes; for the values of t for which uq + Ivq does not belong to K, a Taylor 
expansion shows that 

PdK{uo + tvo) = uo + tvo + 0(t^), 

hence (3J^). With the change of variable U'^ — U^ A-tV'^, equation ( 1.16| ) is equivalent 
to 

V = G{v) 



where the function G is defined by 



G{v) = -V"" + 



1 + e 



Pk U"^ 



2h 

TT"e 



y° 



F[h,U\U 



1 .,0 y°+^ 



M]-u' 



l + e 



Let us check that G is a strict contraction on ^^^(0,2 \vo\^^ + l): if |w|^^ < 2 \vo\^ 
then 



w + V"" 



< 



1 



1 



\Muo + o\^0 + -^ih)\uo + o 



therefore, for h small enough, \{v + Vb)/2|^ is at most equal to 2 |wo|„ , and we can 
use the definitions of L and Ci : 



F{t,,U\U',^^,h 



<Ci+L(2|«oL + l). 



(3.2) 



We estimate G{v) as follows: by the triangle inequality, and the Lipschitz condition 
on Pk, 



|G(.)L„ < i±^7 



U 



2h{V''-vo) 



-y" 



l + e 



Pk 



Uq 



l + e 
2hv{) 



l + e 



l + e 



F-uo 
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We apply (|3j), ( |1.14D and (|3^, and we find 



|G(«)|„„ < 7 2 |z(/.)L^ + ;i(Ci + L{2 |«ol„„ + 1)) 



\vo 



Kh)\ 



"0 2^ 



4/i 



2/i 



Therefore, for h small enough, G maps i?„o (O, 2 |uol„g + l) to itself. Morover, the 
Lipschitz constant of G on this ball is at most equal to 'yLh/2. This proves that G 
has a fixed point in i3„o(0, 2 |uo|„ + l) for small enough values of h and completes 
the proof of the lemma. D 

Here is the statement of the uniformizable estimate, which will be used throughout 
the induction. We will say that t/° and U^ satisfy condition E(il, tq, C2, /i) if 



|t/°-u| <ro, 



\U^~u\ <ro. 



\U^-U°\ <C2h 



(3.3) 



and moreover U^ is uniquely defined in _B(u,ro) by 



1 



and the following inequalities are satisfied: 



|C/^ -u\ <r-o. 



iL/^-L/H <C2/i. 



(3.4) 



(3.5) 



Lemma 3.2. For all u G K , there exists tq such that for all C2 > it is possible 
to find hi > and C'2 < 00 with the following properties : fo r a ll h < hi and for all 
choice of U^ , U^ satisfying co nditi on E(u, rp, C2, h) i.e. (3^), (3^) and (3^), there 
exists a unique U^ satisfying ( 1.16 ) for m — 2 and the estimate 



U' 



U^\ <C'.h. 



Proof. Subtract ( [l.lq ) for m = 1 from ( |1.16D for m — 2; with the change of 
variable U^ — U'^ + hV^ , we have to show the existence of a solution of 



y2 = -eV° 

1 



1 



+ e f2U^ + 2hV^ -{l-e){U° + hV^) + h^F^ 

TT^^l iT-e 



e^ / 2[/i - (1 - e)t/° + /i^F^ 
Pk 



1 + e 

If we denote by G(V'^) the right hand side of the above equation, we have to choose a 
parameter C2 such that G will be a strict contraction of the ball 5(0, C2) into itself. 
Let ro be such that Pr- is Lipschitz continuous on B(u, 2ro); denote by 7 the Lipschitz 
constant of Pk over this ball, and define 



C'2 - (37 + 1)C2; 

let Ci be given by 

Ci =sup{\F{t,u,u',0,h)\ -.t e [0,T], \u~u\ < 2ro, 
\u' -u\ <2ro, he [0,h*]}. 



(3.6) 



(3.7) 
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Denote finally by L the Lipschitz constant of F defined as follows: 

L = sup<^ j : t e 0,T , u-u < 2ro, 

1^ \v — v'\ 

\u' -u\ < 2ro,\v\ <C2, |w'| <C2, he [0,h*], v i^ v' 

Then, we have the estimate for \v\ < C2: 

\F(t2, U\ C/2, {v + V^)/2, h)\<Ci+ LC'2, 
\F{ti, t/°, U\ (V^ + y°)/2, h)\<Ci+ LC'^. 

It is straightforward that 

2L/1 - (1 - e)[/" + K^F^ 



(3.8) 



1 + e 
2[/2 - (1 - e)C/i + /i2i?2 



<ro + 



1 + e 
2hC'2 h^jCi+LC'^) 



1 + e 
Therefore, if hi satisfies the estimate 

2/iiC^ hliCi+LC'^ 



1 + e 



1 + e 



< ''o, 



(3.9) 



we may use the Lipschitz continuity of Pk on the ball of radius 2ro about u, and we 
find that if v belongs to -6(0, C2), 

\G{v)\ < eC2 + 7((3 - e)C2 + 2h{Ci + LC"^)); 

We observe that 7 is at least equal to 1, since DPqk has eigenvectors relative to the 
eigenvalue 1; therefore 

e+(3-e)7 <37+l; 

thus, if h is so small that 

[e + (3 - e)7]C2 + 27/11(^1 + LC'^) < (37 + 1)6*2 = C^, (3.10) 

G maps i?(0, C'2) into itself; moreover, the Lipschitz constant of G over this ball is at 
most equal to ^Lh/2] if 



-jLhi < 2, 



(3.11) 



G is a strict contraction from i3(0, Cj) to itself, which proves the lemma. D 

When u belongs to dK, we need local coordinates in which the projection Pk is 
particularly simple. They are defined in the following fashion: we choose a coordinate 
frame in R'^ such that 

• u — 0; 

• the tangent hyperplane to dK at is the hypcrplane of the first d — 1 coor- 
dinates; 
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• the positive direction of the d-th coordinate axis points inside K. 
For a d-dimensional vector x, we will use the notation 



{xi 



,Xd-l). 



Locally, dK is a graph over the hyperplane of the first d — 1 coordinates, and it 
can be parameterized as 



X{x') 



X 

H{x') 



where x' belongs to R'^~^, H is of class C^ and DH{0) vanishes. Let s i-^ ip{s, z) be 
the parameterization of the geodesic starting at z e dK with an initial velocity equal 
to —N{z) which satisfies 



9V 
9^ 



(s,z) 



i>{s,z) 



Let ^ be defined by 



*(x',2;)=^(-y,x(x')); 



(3.12) 



(3.13) 



the function ^ is of class C^ in a neighborhood of 0; its derivative at has the block 
representation 



L'*(0,0) = 







iV(0) 



(3.14) 



it is invertible, since A^(0) docs not belong to the tangent plane at to dK . Thus 
^ is a local diffeomorphism from a neighborhood W of to a neighborhood 'i!{U) of 
0. In particular, we may assume that U contains a compact neighborhood of of the 
form O X [— ri,ri] where O is an open neighborhood of in M''"^. 

The inverse diffeomorphism of $ is denoted by $, and we decompose it as 

*(x) = (^g) , (3.15) 

V^'^ and Y takes its values in R. If x belongs to 
^=*(Ox [-ri,ri]). 



where 5* takes its values in 



the projection Pk{x) is given by 

Pk{x) = * 



Six) 
Y{x)- 



(3.16) 



With these preparations, we are able to prove the main local estimates: 

Theorem 3.3. For all u G K , for all C2 > 0, there exist two positive numbers, 
ri < r2 and three numbers t > 0, hi > and C3 < 00 such that for all h G (0, hi] 
and all to G [0,T), for all U^ and U^ , satisfying the condition E[u,ri,C2, h) , [Z™ is 
defined in B{u,r2), for all m < [t//iJ, and \V"^\ is bounded by C3 independently of 
h forO <m< [t/H] - 1. 

Proof. The theorem decomposes into an easy and a difficult part. The easy part 
is when u belongs to the interior of K. 
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First case: u G mt{K). We choose rg > as in tlie proof of lemma 3.2: the 
ball of center u and radius 2ro is included in K. The number C2 defined by (3^9) is 
equal to 4C2, and the numbers Ci and L are given respectively by (3/7) and (|3.^). 
We choose ri = ro/2 and r2 = rg. Assume that r satisfies the following inequalities: 



Then, if we write 



< r(Ci + LC0 < min(C^ - C2, tq - n), 
tC2 + ^^{C,+LC'2)<^. 



n = lT/h\ , 



(3.17) 



(3.18) 



we shall prove by induction that for small enough h, there exists a unique solution 
of ( 1.16 ), for < m < n, which satisfies the estimate 

We claim that for h small enough, we can find a solution of 



U^ 



771+1 



2W 



U'' 



= h^F (i„, C/"-\ C/™, {V" + y"-i)/2, h) 



which satisfies the estimates 



Vme {0,...,n-l}, \U"' -u\ < n + Tnh{Ci + LC'^) , 
Vme{0,...,n}, \V"'\ <C2+mh{Ci+LC!^). 



(3.19) 



(3.20) 
(3.21) 



In this construction, we seek a solution without considering the constraints, and we 
prove eventually that they are satisfied. 

It is clear that (3.20) holds for to < 2 and that (B.21) holds for m < 1. Assume 
that it holds up to some exponent m < n. Thanks to (3.17), we have the estimates 



|C/"-2_y| <^g^ 



\U' 



rn — 1 -77 



-^ <r-o, |(7'"-u|<ro, 



\V 



m~2\ 



<c' 



\V" 



<C!,. 



Therefore, we may apply lemma 3.2 with K — M'*: defining C2 — 46*2, we can find 
hi such that ioi < h < hi, there exists a unique [/'"+i such that 



U' 



m+l 



U' 



< C2h. 



In particular, if L' is defined by (3.8), with C2 replaced by C2 , we infer from ( 3.19| ) 
that 

|^m| ^ |ym-i| ^i^^Ci+L'C!^); (3.22) 

therefore, with the help of the induction assumption, we have the estimate: 

|y"| < C2 + mh{Ci + LC'2) + h{L'C'^ ~ LC'2). 

If h satisfies the inequality 

C2 + t(Ci + LC^) + h{L'ai - LC'2) < C'2, 
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we can see that in fact 

|y'"| < c'^, 

and therefore, instead of ( 3.22| ), we have 

IT/"! <C2+mh{Ci+LC'^). 
Therefore, we have also 

Thus, (|^ and ( ^^ hold. Let us prove that the vector VF™ defined by ( jTl^ ) 
belongs to K: since 

Oh h'^ F"^ 

1 + e 1 + e 

we have the estimate 

1 + e 1 + e ^ 
thus, ii h < hi and /ii satisfies 

2/iiC^ 



1 + e 



+ /i2(Ci+LC0<ro, 



H^™ belongs to K, and the sequence f/™ satisfies (1.16). This concludes the proof of 
the estimates in the first case. In particular, we can choose C3 = C2. 

Second case: u e dK. We define on R'' a norm denoted by || || as follows: 



a;= i^] , l|a;|| =max(|a;'| ,\xd\). 



Pick i?i > such that ^ is a diffeomorphism from an open neighborhood 

B={x: ||a;-$(u)|| < i?i} 

to its image and such that 'i'{B) is included in an euclidean ball B{u,rQ) such that 
Pk is Lipschitz continuous on B{u, 2ro); denote by 7 the Lipschitz constant of Pk on 
B{u,2ro). 

Define A by 



r r\D-i'(x)xi\ „ , } 

A =max<^ sup<^ - — r--^. — - : x e B,xi^Q}, 
[ >- Ikill J 



and 



C4 = max<^ sup<^ — -—, — 7--. — : X e B,xij^0,X2 j^Q> 

I '^ 2||a;i|| ||a;2|| J 

I 2\\ui\\ U2 J 
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A continuity argument shows taht the compact set "^(B) contains the baU of radius 

Ri/A about u. 

We win give now a description of the scheme ( 1.16| ), ( 1.17 ) in the new coordinates 
X™ = $([/™). Assume therefore that 



C/' 



Tim I rpTjin — l 
ra+l jjra jj^-l ^m ^^j 1 ^^j ^^ ^ 

1 + e ^ 



We know that (1.16) is equivalent to 

iT'e 



= Pk{W"'); 



(3.23) 



(3.24) 



We map (3.24) by $, and we calcu late the Taylor expansion of either side of (3.24) 
around L/™. The left hand side of (3.24) can be rewritten as 



U" 



' 1 + e 



and therefore 



$ n/™ + /i 



V"^ -eV 
1 + e 



m — 1 



$([/") + D^u"')h + r 



where 



I/" II < C. 



/j(y™_ey™-i) 



1 + e 



But 



SO that another Taylor expansion gives 

$(C/™+i) = $([/") + D$(C/'")/iF" + J™ 
with 

/" <C4\hV"'f. 

Thus 

D<P{U"')hV'" = $(C/™+i) - $([/") - /'". 
A similar calculation gives 

-D'^{u"')hv"'-^ = $(L/'"-i) - $(^7") - r 

with 

<c4l/lT^'"~M^ 



(3.25) 



(3.26) 

(3.27) 
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<^{U 



m+l\ _ 



1 + e 



T 

1 + e 



where 



T =/™ + er"-(l + e)/™, 



and have the estimate 



L,2/'n^m|2 I „ lT^m-l|2 _|_ /^ _|_ g\-l lym _ g-^^m-l|2-. 



<Cih^{\V"'\^+e\V'' 



(3.28) 



Consider now the right hand side of (3.24). By definition of V"^ ^, we have the 
identity 



and a Taylor expansion gives 



(1 - e)hV"'-^ + h^F 



*(W^'") = $([7™) + !?$([/'") 



with 



|J"II <Ca 



(1 - e)hV"'-^ + h^F 
TT'e 



J" 



(1 - e)hV'^-^ + h^F 



1 + e 



We substitute (3.27) into ( 3.30| ), and we obtain 

2X" - (1 - e)X™-i + h^D'^{U"^)F'^' 



$(VF™) = 



J 



1 + e 



1 + e' 



where 



so that 



J^-(l + e)J" + (l-e)r 



<Ca 



We have to estimate 



\{l-e)hV"'-^ + h^F' 



1 + e 



+ {l-e)h'\V 



2 I T/TTT, — 1 I 



|J™|| <C4h^ 



; by elementary inequalities, 



2(l-e)2|y'»~i|%2/i2|F™|^ 



1,2 9 , 1,2 2|V^"f + 2e2|l/™-i| 

+ (l-e) l/"-H +|y"|Ve y"-H +-^ 



1 + e 



The coefficient of | V"^ ^ | in the above bracket is 

2(1 -e)2 2e2 



1 + e 



1 + e' 



(3.29) 



(3.30) 
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and since for e G [0, 1], (1 — e)^ < 1 — e-^, this coefficient is at most equal to 3. The 
coefficient of jT^™! in the same bracket is at most equal to 



1 + e' 



which is also at most equal to 3. Therefore 



< C4/i^[3 |y"r + 3 |y'"-Y + 2/i^ 1^"!] 



Thanks to the properties of P^, 



1 + e ' 






Define 

s = S{U )^ [X \ , y = Y{U ) = A^ 
In these new coordinates, we have 



y^n+1 ^ ^yrn-1 _ ^^^m _ ^^ _ e)y"-l)+ = h' X^ ^ 

and K™ and A™ are given by 

/i^A" = [2y"' - (1 - e)2/'"^i + {h^ D^{U'^)F'^ +7^) ^'^ 
-(22/™-(l-e)y'"-i)++/r. 

Therefore, we have the estimates: 

max(|K™| , |A™|) < A |F'"| + 6*4(3 |y"|^ + 3 jy^-if + 2h'^ |F™|^) 
We define ^™ and C™ by 



C^ 



Let now q be a number which satisfies 



Let 



q > AC2. 



C^=Aq(37+1), 



(3.31) 



(3.32) 



(3.33) 
(3.34) 



(3.35) 



(3.36) 



and let Ci and L be respectively as in (3/7) and (3^). If we assume beyond ( 3.23|) 
that 



max( 1^ 



Tn — 1 1 IT/Tn 



F™i) < a, 



(3.37) 
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we have the estimate 

|F™| < Ci +L(|y'"| + |y™-i|)/2; 

by elementary inequaUties, 

IF"!^ <2Cf + L^{\V'' 
and therefore, if we define 



"2/11^^^12 i_ iT^m— 11 



Cs = ^ + Cl (A + AhjCi) , C-e = ((I + 2hlC4 ) L' + 3C4 ) A^ 



we have shown that under assumptions (3.23) and (3.37), the following inequality 
holds: 



Iini<c5 + c6(iirf + ||r-i ). 

Let T be a number which satisfies the following inequalities: 
r > 0, AC2 + (2rC5 + 2C6q^) < q, 
< p = ^ - TAC2 - 2t2 (C5 + 2Ceq^). 

Assume that initially 

max \W -u\ < Ri/(2A), max \W+^ - W\ < C2/1. 

i=0,l,2 I I j=0,l ' ' 

We will prove by induction that if n = [t//iJ , then for all to < n 

v/e {o,...,to}, ||x'-$(u)|| < ||a:°-$(u)|| 

+ lhAC2 + 2l{l - l)/i2 (C5 + 2C6q^) , 
V/e {0,...,TO-1}, He'll < AC2 + 2lh{C5 + 2Ceq^) . 



(3.38) 



(3.39) 



(3.40) 



(3.41) 



For m < 2, assumptions (3.4C) guarantee that (3.41) holds. The induction hypotheses 
imply that 

max \\X^ ~ $(u)|| <Ri- p, max \\X^'^^ - XH\ < qh, 

j—7n—2.7n—1.7n j— m — 2,m — 1 



and therefore, C/™ ^, [/™ ^ and C/™ belong to ^(S). We may apply lemma (3.2) 
which guarantees the existence of [/™+^ such that 



|C/™+i-C/"| < (37 + l)Ag/i. 
The ball of radius p/A about C/™ is included in *(S); thus, if 

(37 + l)A^qh < p, 
C/"+i also belongs to *(S). Similarly, 



(3.42) 



1 + e 



-[/" 



< (37 + l)Ag/i, 
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(37 + lA^qh < p, 



(C/™+i + eC/™"i)/(l + e) belongs to ^{B). Finally, thanks to the definition ( |3.36| ) of 
C2, we have the inequality: 



U' 



rn+1 



C/" 



2h 



<(^I±^<a 



2i 



in virtue of the definitions ( |3.8| ) of L and ( p.TI ) of Ci , we have 
Once again, if 



1 + e 1 + e 



1 + e 1 + e ^ ' A 

VF™ belongs to ^(6). Thus ( 3.23| ) holds and we may apply the argument that followed. 
By definition of cr™ and 77™, we have the inequalities 



and thanks to lemma |2.1| 



+ h\\C 



^m — 1 1 



■^•Tn — 1 I 



|?7"| <max(|?7'"-^| ,e|7?'"-2|) +/i||C"i +/i||C" 
hence we infer that 

lldl < niax(||C"l >e |I^™"1) + '^ IK"li + ^ l|C""i ' 
and thanks to the induction hypothesis 

IICII < AC2 + 2mh{C5 + 2Ceq^) - hCeq^ 
+ h{C5 + Ceq^) + hCe\\rf- 
The equation in a 

hCea^ -a + 2mh{C5 + 2CQq^) - hCeq^ + h{C5 + C&q^) + AC2 = 
has two distinct real roots if 

A = 1 - AhCQ{2mh{C5 + 2C6q^) - hC^q^ + h{C5 + C^q^) + AC2) 
is strictly positive; but this is always true if < ft- < /ii and 
1 > 4/iiC6(2t(C5 + 2^69^) + AC2). 



(3.43) 



(3.44) 



(3.45) 



Th e sm allest of the two roots of ( p.45| ) is inferior to q, since the substitution a = q 
in (|3.45| ) gives a negative l eft ha nd side; the largest of these two roots is at least equal 
to l/(2/iC6); but relation ( 3.42 ) implies 

liril<AC^; 
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thus, if 



hi < 



1 



2C6Aa 



relation ( 3.44 ) implies ||^'"|| < q: if we substitute t his in equality in the right hand 
side of (3.44), we find that the second inequality in (3.41) holds for I — m; the first 
inequality in (3.41) for I — m -\- 1 holds immediately, and the induction is proved. 
Thus, we can take as an upper bound of |y"| the number C3 = Aq; we can take also 
n = i?i/(2A) and rz = ARi. D 



If we put together theorem 3.3 and lemma 3.1, we obtain an existence result: 
For all {tQ,UQ,M{uo)vo) G 



Theorem 3.4. For all {to,UQ,M{uo)vo) e ED, for all U^ satisfying (1.14), there 
exists T > 0, C3 < 00 and hi such that for all h e (0, /ii], there exists a unique 
solution of (1.16|) and (1.17[) for all m < [r/h] — 1, which satisfies the estimate 



\JKn~l, 



\V'\ <C.. 



(3.46) 



Proof. Let us check that [/° and [/^ satisfy condition E{uo,ri, C2, ft-). Lemma 3.1 
and assumption (1.14) on U^ imply that 



\U'-uo\<h{\zih)l^ + \vo\J 



and 



\U'~uol^<h{2\vol^ + l). 

Choose C2 > (4|wo|„Q + l) ||M(wo)||; U° and U^ satisfy condition E{uo,ri,C2,h) for 
small enough values of h. Then, it is clear that theorem 3.3 applies. D 



It is convenient to give a uniformized version of theorem 3.3: 
Theorem 3.5. For all compact subset C of K, for all C2 > 0, there exist 
positive numbers ri, r2 > ri, t, C3, and hi such that for all to G [0, T), for all 
u G C, for al l h < hi an d for all [/° and U^ satisfying co nditio n E{u,ri,C2,h) re- 
lations ( 1.16| ) and (1.17) define uniquely under condition ( 3.46| ) the vectors t/™ for 
2<m< [min(T,T - to)/h\]. 

Proof. Any element w of C is included in an open ball int(i?(w, ri{u))) such that 
theorem |3.3| holds. We cover C by a finite number of balls mt{B{uj, ri{uj)/2)) with 
associated numbers r2(uj), t{uj), hi{uj) and Cj,{uj). If we let 



ri = - min{ri(uj) : I < j < J}] 

then any M G C belongs to a ball B{uj, ri(uj)/2), and in particular, B(ll, ri) is included 
in B{uj, ri{uj)). If we take 



T = minT(Mj), r2 = maxr2{uj) hi = imnhi^Uj), C3 
j J J 



max C3 {u 

j 



3J' 



it is immediate that the theorem holds, thanks to theorem 3.3. D 



4. Estimates on the acceleration. In this section and the three following ones, 
we assume that there exist strictly positive numbers r, C3 and hi, and a subsequence 
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of times steps to which correspond solutions of the numerical scheme defined by ( |l.l3 ) , 
( 1.14 ), ( 1.16 ) and ( 1.17 ), which satisfy the estimate, for all h < hi: 

V;e{0,P-l}, \U^+^ -U^\<C3h (4.1) 

where 

P = \rlh\ 

Here we estimate the discrete total variation of the sequence (y™) . It is also 
convenient to define the function Wh(t) on [Lq, io + t] by 

Wh{tm) = W™" , Wfi is continuous and it is affine 

on each interval [tm, im+i), and constant on [tp, to + t]. 



Theorem 4.1. Under assumption ( [4.l[ ), there exists a constant C7 such that for 
all h < hi: 

p-i 

Y^ |ym_ym-l| <(^^_ (4 2) 

Proo f. Let C be the compact set K n B{uo^Cj,t) and let ri be as in theo- 
rem ^.5| cover C with a finite number of balls B{uj,ri/4); observe that, thanks to 
Ascoli-Arzela's theorem, the set W of functions {wh)o<hi<h is relatively compact in 
C'^{[to,to + 'T]). The set of limit points of {wh)o<ch<hi as h tends to is also a compact 
set, which we shall denote by Woo- There exists a finite subset w-^, . . . ,w^ of Woo 
such that 

Vw e Woo : infill w ~ w'\\co[t„,to+r] :!<*</}< n/A. 

For each i S {1, . . . , /}, it is possible to find a finite increasing sequence of times 

= T{i, 0) < • • • < T{i, k) < ■ ■■ < T{i, K{i)) = r 

such that 

w\[T{i, k), T{i, k + 1)]) C B(Mj(,,fc), ri/4). 

Thus, for all w G Woo, 

w([r(i, k), r(i, k + 1)]) C B(Mj(,,fc), ri/2). 

Therefore, we can decrease hi so that 

V/iG(0,/ii], 3ie{i,...,l}, Vfce{i,...,K(i)}, 

Vte [T{i,k),T{i,k + l)] Whit) eS(uj(,,fc),3ri/4), 



and thanks to (3.29) and to ( |4.lD , we can decrease hi such that 
V/iG(0,/ii], 3ie{l,...,/}, Vfce{l,...,K(i)-l}, 

yie{lTii,k)/h\,...lT{t,k + i)/h\}, C/'ei?(H,(,,fc),ri). 

We simplify the notations by letting 

P^[T{i,k)/h\, Q=[T{i,k+l)/h\, 

and we take Ci be as in ( p. 7] ), where u is set equal to Uj(^i^k)j i^o is set equal to ri and 
C2 is set equal to C3. 

Now, we have to consider two cases: 
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First case: B(u,ri^ D dK = 0. We have the inequahty 

|F"| <Ci+LC3, 



(4.3) 



hence, thanks to (1.16), we have the inequality 



|^m_ym-l| <h{Ci+LC3) 



and therefore 



J2 |y"-y"-i| < {T{i,k + l) + 2h-T{i,k)){Ci+LC3 



(4.4) 



m=P+l 



Case 2: B{u,ri) n dK ^ 0. We observe that thanks to ( 3.38 ), we have the 
estimate 



V7ne{P+l,...,Q-l}, maxd^^hlA"!) < Cg 



(4.5) 



where 



Cg — C5 + zCg A G3 . 



The estimates on the first d — 1 components of the velocity in the straightened 
coordinates are immediate: 



E 



h 



<{T{i,k + l) + 2h-T{i,k))Cg. (4.6) 



m=P+l 

In order to estimate the last coordinate, we partition {P + 1, . . . , Q} as follows: 
P={me{P+l,...,Q}: 2y™ - (1 - e)y"-i < 0}, 

r'^{p + i,...,Q}\r. 

We write V as an union of discrete intervals: 

£ 

v = \J{p{i),...,q{i)}, Pii)-Hr, q(i) + iir. 
1=1 



If 7/™ is defined as in (2.2), we observe that for m £ T", 

I ?7i _m— 11 ^ i^-f J 
|?7 - ?7 I < Cg/l, 

SO that 

J2 If?" -77"'"'! </iC9|7''|. 

If m belongs to P, we observe that 

??'" - ry™-i = h\"' + (2y'" - (1 - e)jy™-i)" 
and therefore, by the triangle inequality, 

1^™ _ ^m-i| < /^^g ^ (2ym - (1 - e)y™-i)" 



(4.7) 
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and using (4/7) again, 



\V ~ V < ZnCg + Tj —J] 



(4.8) 



Wc observe that we have the elements of a telescoping sum: we sum (4.?) for m 
varying from p{l) + 1 to q{l), and we obtain 



«(0 

Y^ |^'»_ry"^i| <2/iC9(g(/)-p(0)+r?''*'^-r;P('^ 

n=p{l) 



(4.9) 



Now, we sum (4.9) from / = 1 to £, which yields 



I q{i) 

^ — 1 m—p{l) 



< 2Cg{Q -P)- 7?P(i) + r;«(^) - E '/"^'^ " '7^^'"'^ 



;=2 



But the terms 77^*^'^ — 77''^' -'^^ can be estimated, since they correspond to a summation 
over V': 



j^P{l) _ ^9('-l) 



<C9/i(p(0-9(0)- 



Therefore, we have proved that 



E h™-^™"^! <3C9/i(Q-P) + 2C3A. 



=p+i 



Summarizing this relation with (4.6), we can see that 



Y^ |l/™_l/™-i| < AC9(4r(i,/c + l)-4r(i,/c) + 8/i)+2C3A. (4.10) 

m=P+l 



Relations (4.4) and (4.10) do not depend on h < hi; since we have only a finite number 
of these estimates, the theorem is proved. D 

5. Variational properties of the limit of the numerical scheme. In this 



section, we work under the assumption (4.1). Recall that n — [t//iJ. We define a 
function u^ by afRne interpolation, as follows: 



jj-m+l Trrn 

Uh{t) = C/™ + (t- mh)^ ^p^ 

for t G [mh^ (m + l)h), < m < n — 1, 
^Uh{t) = U'' iort(E[nh,T]. 



(5.1) 



We also define a measure Fh as the following sum of Dirac masses: 

71 

Fh{t) = E hF'^Sit - mh). (5.2) 
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In this section we prove that the sequence {u}i)ii converges in an appropriate sense 
to a fun ction u which satisfies (L3) to ( f .5f| ) with r instead of T. We delay the proof 
of (1.6), the transmission condition at impacts, to a later section. 

There are three steps in the convergence proof: the first is to prove that the 
limit u exists in an appropriate sense and takes its values in K; in the second step, 
we show that Uh is of bounded variation uniformly in h and that Fh converges to 
M{u)^^f{-,u,M{u)u) weakly in the space of R'^-valued measures. The last step is 
the characterization of the measure /i = M{u)u — /(•, u, M{u)u): there we show that 
^, satisfies conditions ( 1.4a ), ( 1.4b| ) and ( 1.4cD . 

Lemma 5.1. From all sequence of functions {uh)h indexed by a sequence h tending 
to 0, it is possible to extract a subsequence, still denoted by (m/J^ such that 

(5.3) 

(5.4) 

The function u takes its values in K . 



Uh - 


-^ u in C°([to, to + t]) strong. 


Uh - 


-^ u in L°"{[tfi, to + t]) weak * 



Proof. Thanks to assumption (4T), we know that (uh)o<h<hi is uniformly Lip- 
schitz continuous over [to, to + t]. Therefore, we may extract a subsequence, still 



denoted by Uh, such that (5_^) and ( p.4[ ) hold. Thus u belongs to W 



C'^{[to, to + t]), which means that u is a Lipschitz continuous function 
belonging to {1, . . . , n}, we have: 



to,ta 



in 



For all 



Z" 



u 



711 -\-l 



eU 



m — 1 



1 



[/" 



h 



v- 



■eV" 



1 + e 



(5.5) 



hence C/™ = Z™ — h{V"^ — eV"'' ^)/(l + e). By definition of the scheme, we have 
Z"' = Pk{W"'){(^^), and thus Z™ belongs to K. It follows that, for all m e 
{I, . . . ,n}, the euclidean distance between [7™ and K can be estimated as follows: 



nn{|[/™-M| -.ueK} </i|F"-ey""i| /(1 + e) < /1C3. 



(5.6) 



Thanks to the definition (5.1), we can see that for all t e [io,io + t] the euclidean 
distance between Uh{t) and K is estimated by 2/1C3. This allows us to pass to the 
limit when h tends to and to conclude. D 

Next lemma describes the convergence of the measures involved in our problem; 
we denote by M^((io, ^o + t)) the space of bounded measures over {to, to + t) with 
values in W^. 

Lemma 5.2. The measures ilh and Fh converge weakly in M^({to,to + t)) re- 
spectively to il and M(u)~^ f{-,u,M[u)u). 

Proof. The measure ilh is a sum of Dirac measures on {to, to + t), more precisely, 
we have: 



n 

il^{t) = Y^ {V^ - V'^'^)d{t - mh) - V''5{t - {n + \)h), 

m—l 

and the total variation of iih on {to, to + t) is estimated by 






TV {ilh) < 1^ |V^" 
?n— 1 



(5.7) 



(5.8) 



Theorem LI implies that {uh)o<h<hi is a bounded family in BV {{to,to + t)), the 
space of functions of bounded variation over {to, to + t), with values in R'^. Using 
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Helly's theorem, we can extract another subsequence {uh), which converges, except 
perhaps on a countable set of points, to a function of bounded variation. Hence 

u^BV{{toM + T)). 

Moreover, 

iifi ^' u weakly in AI ((to, ^o + t)). 

Lebesgue's theorem implies that iih converges to ii in L^(to, to + t). We extend 
iih and ii to M by outside of (to, to + t) and still denote the respective extensions by 
u/i and ii. The set {iih ■ h G {0, hi]} U {u} is a compact subset of L^(M). The classical 
characterization of compact subsets of L^(]R) M implies that 



lim sup \uh{t-6)-Uh{t)\ dt^Q. (5.9) 

''^Oo</i<hiJK 

Letting 6 — h, we can see that u/i(- — h) converges to u in L^(M). Let us define an 
approximate velocity v^ on R by 

The sequence Vh converges to ii in i^(lR). Moreover, for all t e [tm,t„j+i) and for all 
m G {1, . . . , n}, we have the identity 

yTTi , ym-l 

vh{t) = -^ . (5.11) 

We have immediately the following estimates for all t g (to, to + r) and all h S (0, hi]: 
\vh{t)\ < C3; \uh{t) - uol < Csit - to) < C3T. (5.12) 

Let V' be a continuous function over [0, T] with compact support included in (to, to+r). 
For all small enough h, the support of ip is included in [to + /i, to + nh]. The duality 
product {Ffi , f/;) has the expression 

n 

{Fh, V') = ^ hi;{to + mhfF"'. (5.13) 



We wish to compare the expression ( 5.13 ) to 

rto+T 



ilFM{u)-^f{-,u,M{u)u)dt. (5.14) 

tn 



We compare the right hand side of (5.13) which is basically a numerical quadrature 
by the formula of rectangles to an appropriate integral. Let us rewrite the individual 
terms of the right hand side of ( 5.13| ) as 






hi;{tmfF"'^ i;{tfF'''dt+ {^{t^) - ^{t)f F^ dt. (5.15) 
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Consider the second term on the right hand side of ( p.lq ): we have aheady proved 
(see relations (13) and ( |4.5| )) that there exists a constant Cg independent of m and 
h < hi such that 



max |F" 

0<?7i<n 



<Ch 



Denoting by w^ the modulus of continuity of ip we can see that 



{^{t„,) - ^{t)f F^ dt 



< CsUJ^{h)h. 



(5.16) 



(5.17) 



We consider now the first term on the right hand side of (5.15), which we would like 
to compare to expression ( ^.14 ). Thanks to the consistance assumption (1.12) have 
the following inequalities, for all t e [tm, tm+i), and all n G {1, . . . , n}: 

\F'^ - M{uhit))-^ l{t,un{t),M{uh{t))vn{t))\ 

< \F{Un,U"\U"'-\vhitm),h) ~ F{t^,Uh{t),Uh{t),Vh{Un),h)\ 
+ \F{tm, Uh{t), Uh{t), Vh{t,n), h) - F(tm, Uh{t), Uh{t), Vh{t),0) I 
+ \M{uh{t))-^[f{t,n,Uh{t),M{uh{t))vh{t)) - 

f{t,u,,{t),M{u,M)vhit))]\. 
For all t E [to , to + t] , let us define 

Phit)=M{uh{t))vh{t). 
Denote by V the set 

V = |(t, ui, U2, v,h) : < t < T, |ui — iio| < C3T, 

|U2-U0|<C3T, |«|<C3, 0<h<hi}. 

Let L be the Lipschitz constant of (ui, U2) •— > F{t, ui,U2, v, h) restricted to 2? and let 
Up be the modulus of continuity of _F on P With these notations, we can see that 



IF"' - M{uh{t))-^f{t,Uh{t),M{u^{t))u^,{t))\ 

< L(|C/" - uhit)\ + |C/"-i - Uh{t)\) + 2ujF{h). 



(5.18) 



Since M is of class C^ in R'', {uh)h converges strongly in C°([to, to+r]) and {vh)h con- 
verges strongly to ii in i^(R) and almost everywhere on (to, to+r), the sequence {ph)h 
also converges strongly in L^(R) and almost everywhere on (to, to + t) to p = M{u)u. 
We see that M{uh)^^ f{-,Uh,Ph) tends to AI{u)^^f{-,u,p) strongly in_L^ (to, to +t) 
and almost everywhere on (to, to + t). We summarize relations (5.17) and (5.18) 
together with the above convergence result, and we find that 



(Fh,^) 



ta+T 



^^M{u)-^f{-, u, M{u)u) dt 



to+T 



< / \M{uh)-'f{-,Uh,M{u,,)vh,)-M{uy'f{-,u,M{u)u)\ j^M dt 

Jto 

+ Csto^{h)T + (3XC2/1 + 2ujF{h)) [ IV'l dt, 

Jo 
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which concludes the proof. D 

Let us prove now that the measure ^ has the required variational properties: 
Lemma 5.3. The measure fi satisfies properties (L4a), (1.41:) and (1.4c). 
Proof. Define 

Hh = M(uh){uh - Fh); 
fih is a sum of Dirac measures on (to, tg + t); more precisely 



Hh^J2 ^■^(t^")(^'" - V""'^ - hF"')S{t - mh) 

771—1 

-M(C/"+i)y"(5(t-(n + l)/i). 

With all the previous results, we know that fih converges to /i = M{u)il — f{-,u,p) 
weakly in M^({to,tQ + t)). Let us prove property (1.4a). Assume that tq is a point 
of {to, to + t) such that u{to) belongs to the interior of K. Then, by continuity of u, 
there exist e > and p > such that 



,i{\u{t) 



\t-To\ <e,x e dK} > 3/9. 



Since the sequence \Uh) , converges uniformly to u as, h tends to 0, we can decrease 
hi so that 

mi{\uh{t) -a;| : |t-To| <e,0 < /i < hi,x G dK] > 2p. 



Relation ( 1.24 ) implies the identity 



VF" = U" 



1-e 



hV" 



.p" 



(5.19) 



Relations (5.19) and (5.16) imply that 



W 



c/™i <ha 



1 



1 



1 + e" 



(5.20) 



Possibly decreasing hi, we have thus 

ini{\W"^ -x\ : |t„-To| < e,0 < h < hi,x € dK} > p. 

This proves that the support of ph does not intersect the open set (tq — e, tq + e), 
and therefore, relation (1.4a) holds. Assume now that ui = u{ti) belongs to dK, and 
let B(ui, i?i) be a ball having the properties of theorem |3.3| ; assume that the image 
of (ti,T2) by Uh and Wh is included in this ball for all small enough h. We rewrite 
conditions ( 1.4b| ) and (1.4c) as follows: for all continuous function ij} with compact 
support included in (toj^o + t) and taking its values in W^ the following implication 
holds: 



yt^{to,to 



d(t){u{t))'ilj{t) > =^ {p,Tp) > 0. 



(5.21) 



In particular, if d(j){u{t))'ip{t) vanishes for all t £ {to , to + r ), th en {p, t p) als o vanishes. 
The reade r wi ll check the equivalence of ( 1.4b| ) and ( 1.4c ) with ( 5.21 ). We infer 
from relation (5.6) that 

|r(C/™)| < A/iCg; 
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the above relation together with ( 5.20| ) imply that there exists a constant Cio such 
that 

\YiW"')\ < hCw 



Since ( 5.21 ) is local, it is enough to check it in the neighborhood of any ii e (^07 to+r). 
Let 



P^\n/h], Q=[T2/h\, 



and 



V = {me{P,...,Q]:W'-iK], V = {P, . . . ,Q]\V. 
We observe that if m belongs to V' ^ then 

ym _ ym-1 _ f^pm ^ Q 

Therefore, we have the identity: 






m=P 



meV 



We recall relation (1.29). Relation (3.32) implies that 

f 



$(Z") - $(VK™) = 



lr(VK" 



and therefore 



Z™ - VK™ - D^{W''') 



( ' 



<Ci\Z"' -W" 



On the other hand, the definition of ^ is such that the d-th column of I?^(Z™) is 
equal to A^(Z"); therefore 

< 2C4|Z"'-VF™|y(H/")-. 
We infer from the above estimates that 

l^m _ 14/™ __ y(pT/™)-jv(Z")| 

< C4{2Y{W"')- + |Z™ - VK™|) |Z™ - W 

< -^ ly™ _ y™-i _ /iF"| (2 + A)hCi 
1 + e ' ' 

and thus, there exists Cn such that for all ni e P: 

i^m _ ^m _ r(yy™)-iv(Z")| < h^Cn |y'" - y"~i - hF'' 



10, 
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We can see now that 



1 + e 



h 



^(Z"_Ty",V(i™))c 



^^rT. Y{wn-{N{z^)Mtm))i 



mev 



-Ciih{l + e) max (||M(C/'")|| IV-fi™)! V |F™ - F™-i - /iF™| , 
^ mev 

which inipHes by a strightforward passage to the hmit that (m, V'} is non negative. 
This concludes the proof of the lemma. D 

6. Transmission of energy during impact. The basic assumption is still the 
one made at the beginning of Section 0. 

Let T € (0, r) be such that u(r) belongs to dK . Write f = fg + r. We decompose 
p(t ± 0) into a normal component pN^ ± 0) belonging to M.d<t){u(t)) and a tangential 
part pT(t±0) belonging to the orthogonal oi d(j){u(i)) in the cotangent metric at u(t). 

In this section, we shall prove that 

PT(t + 0) = pT(t - 0) and pN(t + 0) = -epN(t - 0), (6.1) 

where e is the restitution coefhcicnt of the problem. 

The conservation of the tangential component of the impulsion is proved in next 
lemma: 

Lemma 6.1. Assume that t e (0, r) is such that u(t) belongs to dK. Then 

PT{t + 0) ^PT{t~0). 



Proof. Thanks to lemma |5.3| , we know that 

M{u)u = ^l + fi■,u,p), (6.2) 

and that there exists a nonnegative measure A such that 

^i = Xd(j){u). (6.3) 



We take the measure of the set {t} by the two sides of (6.2), and we find that 



M{u{t)){u{t + 0) - u{t - 0)) = fi{{t}), 

which implies immediately that p{t + 0) — p{t — 0) is parallel to d(j){u{t)) and proves 
the lemma. D 



Let u = u(t) and let B{u,ri) and _B(u, ri) have the properties of theorem 3.3. 
There exists an interval [T_5,r2] containing r in its interior such the for all small 
enough /i, u;i([to + t-5, to + ^2]) is included in i3(wi,ri). 

The apparently strange notations t_5 and T2 have been chosen in view of the 



upcoming construction of lemmas 6.3 and 3.4, where we will consider relative times 



r_5 < • • • < r_i < T < Ti < T2. 
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Define 

P=rr_5//il+l, Q^\T2/h\-l, 
and let Xh be obtained from the X™ by affine interpolation, for P < to < Q. We 



infer from estimates (4.1) and (4.2) the estimates 

vm+l vr. 



max 

P<m<Q 



E 

m=P 



X 



m+1 



h 

-X' 



<AC3, 



h h 

Therefore, we have the following convergences 



< AC7. 



Xh~^ X strongly in C"([io + t_5, io + ^2]); 

Xh ^ X except on a countable set and weakly * 

ini°°([to + r_5,to + T2]); 

Xh ^ X weakly in M^ ([ig + t_5, <o + ■''2]) • 



Write for all /i < hi 



Xh 



where the ShS and s take their values in W^^^ and the y/i's and y are real valued 
functions. We do not have Xh = ^{uh), because Xh is a linear interpolation of the 
sequence X™ = <l>(f/'"), and $(u/i) is the image of the linear interpolation of the 
sequence C/™. However, we can estimate the difference Xh — $(u/i). 

Lemma 6.2. For all t E [to + t_5, io + T2], belonging to [tm, tm+i], we have: 

Xh{t) -^{uh{t)) < 2C4C|/imin(i-i„,i„,+i -i). 



Proof. We observe that 



Xh y^n 



X' 



and that 



j^[xH{r)-Huh{t))] l^^^^^ 

$(t/'"+i) - $({7™) - /ii)$([/™)y'^ 



h 



< held. 



Moreover, for all t e [t„i,t„j+i 
ci2 - 



df^ 



Xh{t) - '^{Uhit)) 
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Therefore, a straightforward integration yields 

\xh{t) - ^{Uhit)) I < ClC4{h{t ~ tm) + {t- t^f), 

which imphes 

\xh{t) - '^{uh{t)) I < 2C^Cih{t - t„). 

We can write the analogous estimate on the interval [t^tm+i), which concludes the 
proof. D 

As a consequence of lemma |6.2| we obtain: 

yte[to + T^5,to + T2], a;(i) = $(u(t)), 

and 



In virtue of relation (3.14), 

We can rewrite this relation in terms of p^r and px'- 

PT{t ± 0) = M(0) ('^* ^ ^A , PNit±0)^ yCt ± 0)M(0)7V(0) 



Lemma |6.1| implies s{t + 0) = s{t — 0). In order to achieve the proof of relation (|6.1|), 
we will prove the scalar relation 

y(t + 0)^-eyit-0). (6.4) 

We will do t his by performing a precise analysis of the transmission of energy on the 
scheme (|2.l|) . The measure jjh is a sum of Dirac measures on (to + t_5, to + T2). We 
define two measures uj^ and A/j on {to + t_5, tp + T2) by 

^ {-2y"^ + {1 - e)y"^-')+ ^^^ 
LOh = > ^ ; — o[t — mh), 



and 



We have 



Q 

Xf,{t)= J2 hX"'Sit-mh). 

m=P 



jjh =UJh + A/i, 



and it is obvious that LUh is a non-negative measure. 

Since the real numbers A™ are bounded independently of h and n, the measure 
by jA^I of any subinterval [a,b] of (to + T^5,tQ + T2) is bounded by C{h — a + h), 
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and it is clear therefore that there exists a function A e L°°{tf) + T_5,to + T2) and a 
subsequence A^ converging to A in the weak topofogy of M^((to + t^5, ta + 'i'2))- 

The measure ujh converges in the weak topology of Af^((io + t_5, io + '''2)) to a 
non- negative measure uj, and in the limit 

y = uj + \ (6.5) 

while 

|A|l=o < Cg. (6.6) 

Since y is non-negative on (ip + t_5, tp + T2) and 2/(to) vanishes, we must have 

2/(Z + 0)>0, y(<-0)<0. 

On the other hand, y(i + 0) — y(t — 0) is equal to uj{{t}); if (^({f}) vanishes, we have 

y(Z + 0)-y(f-0) = 0, 

and the identity 

yit + 0) = -eyit-0) 

holds. Therefore, the only interesting case is when 

cj{{t}) > 0. (6.7) 

The following two lemmas enable us to prove in two steps that the velocity is 
reversed according to the law described by (|l.6|). Lemma 6.3 shows that if uj has 



a Dirac mass at t, then the left velocity at t is outgoing; Lemma 3.4 shows indeed 



that ( |1.6| ) holds. 

Lemma 6.3. If uj{{t}) is strictly positive, then y(t — 0) is strictly negative. 
Proof. The idea of the proof is to find two succesive times im-i l^tm <'t for which 



we can write down an estimate on the discrete velocities, and then to use lemma 2.1 
to perform a discrete integration and to obtain a contradiction. We must deal with 
the fact that yu does not converge uniformly to y. 

Without loss of generality, we may assume that ij is continuous on the right and 
that for all h < hi, i/h is also continuous from the right. According to Helly's theorem, 
there exists a countable set D such that 

yiiit) -^ y{t), Vt such that t-te {t^5,T2) \ D. 

Assume that y{t) vanishes; therefore, y(t + 0) is strictly positive. Choose a = y{t + 
0)/4, and let t_4 and ti be such that 

T-.5 < T_4 < T < Tl < T2 

6C9(Ti-r_4) <a, (6.8) 

and 

Lu{[to+T^i,t)) <a, uj{it,to + Ti])<a. (6.9) 
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An integration of ( |6.5[ ) on appropriate intervals yields 

Vte(to + r_4j), \y{t±0)\ <a + Cg{t-t), 
yt e (t, to + n), y{t ± 0) > Lj{{t}) -a~ Cg{t - t). 



(6.10) 
(6.11) 



Choose T_3 G (t_4, t) \ D and r_i G (t_3, t) \ Z?; since ujh is a nonnegative measure, 
we have the following inequality for all r' G (r_3,T_i) and all t" G (r',T_i): 

IVhito + t') - ijhito + t")\ < Uhiito + t', to + t"]) + Cq{t" -r' + h) 
< iOhiito + T-3, to + r-i]) + C9(r" - r' + h). 

We integrate ujh — w on the interval [tg + t_3, io + ''"-il ; since the measures w and ujh 
do not charge to + ts and to + t-i, we find that 

w/i([to + r_3,to + T_i]) -uj{[to + r_3,to + t_i]) 

= y/i(to + T-i) -yh{to +T-3) -y(to + r_i) +2/(to + r_3) 
+ A([to + r_3, to + t_i]) - A,,([to + r_3, to + r_i]), 

and therefore 

t^/l([to + T-3,to + T_l]) 

< w([to + T_3, to + r_i]) + |y,,(to + T_i) - y(to + t_i)| 
+ \yh{to + r_3) - y{to + r_3)| + C9(2(r_i - t^s) + h) . 

Choose now r_2 G (r_3,T_i) \ D; then, for /i small enough, t^ = h\T2/h\ and 
tm-i = tm — h belong to the interval (t_3, r_i), and therefore. 



\yh{tm) ~yh{tm-i)\ < a + C9(2(r_i -r_3) + 3h) + en, 



(6.12) 



where e^ tends to ash tends to 0. On the other hand, 2)h(to+T_2) tends to y(to+r_2) 
and therefore, thanks to relation ( |6.10 ), there exists a family ej^ such that 

\yhito+T-2)\ ^ |y/i(tm)| < a + Cg{T-T^2) +4> 

which is equivalent to 

It?"! <a + C9(r-T_2)+ek; (6.13) 

we infer from ( 6.12 ) and ( 6.13| ) that 

|77™-i| <2a + C9(2(r_i - t^s) + r - r_2 + 3/i) + e^ + e'^. 
Thus, for all n > m we infer from Lemma |2.l| that 

Iry"! <2a + C9(2(T_i-T_3)+3/i 

+ T - T_2 + 2(t„ - t„j)) + Eh + e'f^. 

Therefore, in the limit, for all t > to + r_2 

\y{t)\ <2a + C9(2(t_i - r_3) + r - r_2 + 2(t - r^a))- 
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and for all t G [to + ^-2,^0 + '''1] 

\y{t)\ <2a + C9(2(r_i - r^g) + r - t_2 + 2(ti - T-2)) ■ 
On the other hand, relation ( |6.11 ) implies that for all t e (t, io + ti), 

\yit)\ >3a-C9(ri-T). 



(6.14) 



(6.15) 



Under assumption (p.q), relation ( |6.15D contradicts relation ( |6.14| ). D 

We can conclude now the local study of the reflexion of the velocity by the fol- 
lowing lemma: 

Lemma 6.4. //aj({t}) is strictly positive, then 



y{t) = -ey{t-0). 



(6.16) 



Proof. Since y{t — 0) is strictly negative, there exists a real number t_3 such that 
y{t) is strictly positive on [to + t_3,Z) C [to + T_5,t). For all r_2 G ('''-a,''^), there 
exists T_i e (t_2,7^) and /ii > such that 



V/ie(0,/ii], Vie [to+T_2,to + r_i),2;,,(t)> 



y(io + T-2) 



(6.17) 



We prove now that there exists a maximal integer 

me{[T_3/h\,...,liTo+e)/h\} 
such that 

yi e { [T_3/h\ , . . . , m - 1}, 2y' - (1 - e)2/'-i > 0, 
and denoting 

Cft = tm-l ~ to, 
the time ah satisfies 

lim ah = T. 



(6.18) 
(6.19) 
(6.20) 



Let us first observe that for all small enough h and all ti belonging to [to+T_2, ^o+'''-i] 
we have 



2y' - (1 - e)y'-^ > 



(6.21) 



Indeed, 



i-i 



2y' - (1 - e)y'-^ - (1 + e)y' + (1 - e)hi 

> ^y(*o + r-2) - /^(l - e)AC3, 



and if 2AC3(1 — e)h < (1 + e)y{to + T-2), we can see that (3.21) holds. Therefore m 
exists and 



lim inf ah > t. 
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On the othe r han d, if there existed ti > r such that for aU t„j £ [to + ^-3,^0 + 
Ti] we had (3.21), then w/j would vanish on (ig + T_3,to + ti), which contradicts 
assumption (3.7). Therefore, we have shown that 

hm sup ail < t, 

i.e. ( |6.20 ). We integrate discretely equation ( |3.33 ), and we find that for t e [to + 

T-3,to + <^h\ 



yh{t) =yh{to + o-/i) - (to + 0-/1 - t)yhitQ + ah) 
Ah((s,to + <Jh])ds. 



I 



(6.22) 



In the limit we have, 



y{t) = y{T)-{T-t)limyhito + c7h + 0)+ / \{r)drds. 



The comparison of ( 6.22 ) and ( 6.23| ) shows that 



limy,, (to + <7h + 0) = limr/'"-i = y{t - 0). 

hlQ hiO 



(6.23) 



(6.24) 



Our purpose now is to obtain very precise estimates on the behavior of y^ beyond 
to + cT/i. Thanks to the maximality of of m, we have the relation 



y 



m+l 



-ey 



771 — 1 I l2 \ m 



+ /i^A"; 



(6.25) 



let us estimate 2y™+^ — (1 — e)y'": we substitute the value of y'"+i given by (6.25) 
into this expression, and we also use (3.34) with m replaced by to — 1; we find 



22/'"+^ - (1 - e)y"' 

= - [2y"'-^ - (1 - e)y'"-2] - (1 - e)h^X"'-^ + 2h'^y 

We apply relation (pTJ) for n = ?7i + 1 and we find that 



(1 - e)/iA™-i + 2hX^y 



m+l _L ct-,'^''^~^ h ( \'^''^^^ X"^~^ 



7/ ' + er/ 

+ {-[2y"-i-(l-e)y 

Therefore, we have 



^-2^J^-l 



77™+! + er;"-i > -2/iC9 



On the other hand, if ^ = -[2y"-i - (1 - e)y"'-'^]h-^ - (1 - e)h\"'-^ + 2/iA™ is 
lesser than or equal to 0, 



|rr+i+e?7™~H <2/iC9 



if ^ is positive, then the sign condition on 2j/™ — (1 — e)j/™ implies that 



Tj^+i + ery"-^ < /i(A™+i + A") - (1 - e)/iA"-^ 
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Thus, we have shown that 

jr^^+i+ery"-!! <3C9/i. (6.26) 

If e is strictly positive, then for all small enough ft,, 

rr+^ >e|y(f-0)|/2. 
Let us estimate now the expression 2y™+^ — (1 — 6)2/™+^: we have 

2ym+2 _ (^ _ e)y"'+^ = -e[2y" - (1 - e)y"'-^] + Oih^). 
If 22/™+^ ^ (1 ^ e)y™ is non-negative, then 

ym+3 ^ 2y"+2 - y™+l + ft2^™+2. 

We must estimate 22/™+^ - (1 - 6)?/™+^: 

2y™+3 _ (1 _ e)y™+2 - 2j/™+2 + (1 - e)y™+i 
= /i(2r/"+2 - (1 - 6)77"+^) 
= ft(l + e)?7"+i + 2ft2A™+2, 

and therefore 2y™+'^ — (1 — e)j/'"+2 is non negative for all small enough /i; the repetition 
of the argument shows that there exists 6* > such that for all small enough h and 
all n £ {ra + 2, . . . ,to + \9/h\}, the expression 2y™'^^ — (1 — e)j/™ is non negative, 
and thus we have the relations 

m— 1 

J/™ = 2/"+i + ft(n - m - 1)77"+^ + ^ (n-.7>2A^'. 
On the other hand, if 2y'^^^~^^ ~ (1 ~ e)y'^^ is negative, we must have 

1 + e 
and therefore 

y--' = -^^ + Oie). 

1 + e 
These relations and the assumption on the sign of 2y™+2 — (1 — e)y"^ imply that 

2ym+3 _ (1 _ ,)y™+2 ^ __(4e^+e(l-e)^)ft<;^ ^ ^ ^^^ (g 27) 

1 + e 

which is strictly positive for h small enough. But now, we can see that 

ym+3 _ ym+2 ^ _ehj^m-l ^ (^(^2)^ 

which is strictly positive for small enough h, and therefore 2y™^* — (1 — e)y™^^ is 
strictly positive for h small enough, since 

2ym+4 _ (^ _ e)y"'+^ > -he{l + e)ry"-i + 0{h^); 
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the same argument as above shows now that there exists 9 > such that for all 
n e {m + 3,...,m+ [9/h\}, 

771—1 

y" = j/™+2 + /i(n - m - 2)77"'+2 + ^ {n-j)h^X^. 

j—ni+3 

If we let a'f^ = tm+i — ^o in the first case and a'f^ = i)ri+2 — to in the second case, we 
have now for cr'f^ < t — to < a'l^ + 6 — h 

yh{t) = yh{to + <yh) + {t-<y'h-io)yh{to + <j'h)+ I A,,((s,t])ds (6.28) 

and 

yh{to + n'^)^0{h), yh{to + n',^)^-eii^'^ + 0{h). (6.29) 



Passing to the limit in ( p.2q ), we can see that 

y{t + Q) = -ey(J-Q). 
If we assume now that e vanishes, relation (|6.26) implies 



7^"+^ = Q{h). 



We observe that lemma 2.1 implies that for all n 

h"| <|V""'|+2C9/i, 
which implies immediately that for n > to + 1 

kri < h'"+^| + 2/iC9(n - TO - 1), 
which proves by a straightforward passage to the limit that 

y(i + 0) = 0. 

This completes the proof of the lemma. D 

7. Initial conditions. In this section we prove that the solution that we have 
constructed satisfies the initial conditions; we work under the hypotheses stated at 
the beginning of section |j. 

Lemma 7.1. The function u satisfies the initial conditions 

m(0) = uo, p(0 + 0)=po- 



Proof. By uniform convergence of Uh to u, it is clear that u(0) is equal to uq. 
There remains to show that the initial condition on the impulsion is satisfied. 

Assume first that uo belongs to the interior of K] then there exist /ii > and 
Ti > such that for all h G (0, hi] and for all t G [0, n] 

\uh{t) - wo| < - inf{|wo ~y\ -y t K}. 
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Then for all i,„ belonging to (0, n], (2C/™ - (1 - e)C/™-i + h'^F"')/{l + e) belongs to 
K for h small enough; we have indeed 



2[/" 



(1 - e)C/™-i + h^F 



1 + e 



uo 



< ^/iC7i + \ inf{|uo - Z/l : 2/ ^ X} + t^Q, 
1+e 2 1+e 

which is strictly inferior to inf{|uo — y\ '■ y ^ K} for h small enough. Thus the 
constraints are not saturated for < tm < ti and the convergence is a classical result. 
In the second case, uq belongs to dK] we have taken admissible initial conditions, 
so that 

(po,d0K)):o>O. 

We use the construction and notations of section 0: $, ^, X™, s™, y™ and C,"^ have 
the same signification as there. 
Taylor's formula yields 

and the definition (1.14) of [/^ gives 

C° = i?$K)Af(uo)-Vo + 0(/i). 
Write 



(7.1) 



D^{uo)M{uo)-^Po. 



Then the normal and tangential components of the impulsion are given by 



Pot = M{uo) ( ° ) and poAf = mM{uo)N{uo). 



We wish to prove 



which is equivalent to 



p(0 + 0)=po, 



i(0 + 0) = 



We recall relation (3.33). Relation (^J) implies that 

a' = {D<i>{uo)M{uo)-'poy + 0{h), 



and together with (3.33), we obtain in the limit 

s{t) = {D<^{uo)M{uo)-'poy + 0(i), 
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i. e. 

s(0 + 0) = ao. 
Let us show now that 

considering two cases: tjq > and tjq — 0. When tjq vanishes, we have 

y'^y" + h,jo + Oih^)^Oih^), 
and 

y' = -ey° + {2y'-il-e)yy+h'X' 
^2{y^f + h^X^^0{h^). 
Thus, 

V" = 0{h), v'^0{h), 
and relation ( ^.3D imphes 

|r/™| <0{h) + 2C9h{n~l); 
therefore, a passage to the hmit gives immediately 

y(o + o) = o. 

If, on the other hand, r]^ is strictly positive, then 

2y^ - (1 - e)2/° = 2y^ = 2hrf + 0{h^) 

which is strictly positive if h is small enough. Let {1, . . . , to} be the maximal interval 
such that 

2y" - (1 - e)y"-i > 0, if n < to. 

Then, for all n G {1, . . . , m], 

which implies by discrete integration that 

??" > Vo - hnCg, 

as long as n belongs to {1, . . . , to}. Moreover, if we choose any ti < rjQ/{2CQ) and if 
n is at most equal to min(TO, [ti//iJ), we can see that 



y™ = yo + /i(^° + • • • + r;^ 



js hnr], 



> 



'0 



2 ' 
for all small enough values of h. 

In particular, for all n < min(TO, [ti//iJ), 

o m n ^ m-l ^ (l + e)/in77o 

22/ -(l-ejy > ^ (l-ej/iACs, 

which proves that to, is at least equal to [ti//iJ . Therefore, LOh vanishes on the interval 
(0, Ti — h); in the limit, lu vanishes on (0, ri) and therefore 

which completes the proof of the lemma. D 
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8. A p ri ori e s tima te s. In t his se cti on w e prove that solutions of the problem 

(Q, ( [l.4a| ), ( |l.4bD ,(|r4^), ( |l.5a| ), ( |l.5bD , ( Q and (U) satisfy an a priori estimate 
on an interval with non empty interior. 

Lemma 8.1. Let R be strictly larger than |p o|,*^. Then there exists t{R ) > such 
that for all solution u of (U), ( |l.4a| ), (|l.4b[ ),([L4^), ( |l.5a| ), ( |l.5b[ ), (^) anrf (|l|) 
defined on [t^, t^ -\' t\, the following estimates hold: 



Vte [fo,to+min(r,T(i?))], \u{t) ~ ua\ < R, \p(t)\l^^■^ < R. 



•1) 



Proof. The measure A appearing in ( 1.4b ) can be decomposed in the sum of an 
atomic part \a and a diffuse part A^. At each point of the support of Xa we have 

b(t+o)i:(,)<b(f-o)i:(,) (8.2) 

thanks to relation (|1.6|). On any interval (^1,^2) which does not intersect the support 



of Aa, we multiply relation (1.3) by v?" on the left, and we find that 
d 1 



Define 



^^;^u^M{u)u = u^f{-,u,p) + -u^{DM{u)u)u. 



E{u,p)^ -(p,p)l, z=\p\l. 



.3) 



It is convenient to recall that 

\p\:^\M{ur'/M^\M{uy/\\. 



Relations { p. 4 ) and (S^) imply that in the sense of measures 

zz ^ E < iF f(-,u,p) + -u^{DM{u)u)u. 



(8.4) 



Our purpose now is to transform (8.4) into a differential inequality. Let xW) be the 
norm of the bilinear mapping 

With this definition, 

\iF{DM{u)u)u\ < x{u)z^. 
We write now 

u^f{t, u,p) = iFM{uf''^M{u)''^/'^f{t, u,p) 

= u^Miuy^^[M{u)-'^'f{t,u,p)-Miuo)-'^^fit,Uo,Po) 
+ M{uo)-'/^f{t,uo,Po)]- 

Define 



g{t) = \M{uo)-'/^f{t,uo,po)\ 
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and let uj{t, R) be the Lipschitz constant of {u,p) ^^ M{u) ^^^f{t, u,p) for t £ [to, to + 
t] and niax(|M — moMpIu) ^ -^i more precisely 



uj{t, R) = sup 



\M{ui)-^/\f{t,Ui,pi)-M{u2) ^'^f{t,U2,P2)\ 



to <t < to + T,max(|ui -uo\,\u2 - uo\,\pi\l^,\p2\lJ < R, 
ui ^ U2 or pi ^ P2 

By construction, uj is continuous and it is an increasing function of r and R. 

Fixi?> |poi:„. 

If to < t < to + T and if inax(|u(t) — mo|, b(^)l^(-t-)) < -R on [to, to + t], we have the 
inequality 

\u^f{-,u,p)\ < z{g + uj{T,R){\u-uo\ + \p-po\))- 

But we can estimate u(t) — uq- 

\u{t)-uo\ < I \u{s)\ds< I \\M{u)-^'^\\zds. 

J to J to 

Therefore we have the estimate 

\u^fi-,u,p)\ < zg + zu;{T,R)( f WMiuy^/^Wz ds 

+ \\M{u)'^^z+\po\y 
and we conclude that z satisfies the differential inequality 

z<9 + uj{t,R)\ f \\M{u)-^^^\\zds+\\M{uy/^\\z+\po 



Xiu)z^. 



Set 



a = sup{||M(M)^/2|| : \u - uo\ < R}, 
(3 = snp{\\M{u)-'^/^\\ : \u - uo\ <R}, 
7 = 2sup{x(m) : \u - uo\ < R} . 



(8.5) 
(8.6) 



While t < to+T and max(|u(t)— uo| , \p{t)\u(t)) 1^ R-t '^ satisfies the following differential 
inequality 



z < g + uo{t,R) (3 I z ds + az + \pt 



- 7Z 



(8.7) 



Let p be any positive number; consider the integrodifferential equation 
y=^9 + p{PJ y ds + ay + \po\j + -i\y\^ , 



with the initial condition 



y{to) = z{to). 
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It has a unique maximal solution which blows up in finite time, as soon as 7 is strictly 
positive and p|pol + sup|.g| is strictly positive. Let 0{p) G [to^T] be the largest time 
for which 

yte[h,e{p)], yit)<R, (3 f yds<R. 

Jto 

As is a decreasing function of p, there exists a unique t{R) such that 

0iLuiTiR),R)) = T{R). 
Choose now 

P = u;{t{R),R). 



Then we can compare the solution z of (8.7) and the solution y of (B.8), and we find 
immediately that 



yt e [to, to + min(T, t{R))] , z(t) < y{t). 
This concludes the proof of the lemma. D 



(8.9) 



9. Global results. We summarize the results obtained so far in the following 
Proposition: 

Proposition 9.1. Assume that there exist strictly positive numbers t, C3 and 
hi > 0, and a sequence of solutions of the numerical scheme defined by (1.13), (1.14), 
(1.16) and (1.17), which satisfies the estimate (4.1). Then it is possible to extract 
from t he seq u ence U h defi ne d by ( |5.ll) a su bseq uence whic h converges to a solution of 
(@, (|l.4a|) , ( |l.4b|) ,( |l.4c[) , (|l.5a|) , (|1.5b| ) , dlj) and (|l|). The convergence holds in 
the following sense: Uh converges uniformly to u^ on [to, to + r]; iih converges to it 
in L°°(to, to + t) weaky star and almost everywhere on [to, to + t], and ilh converges 
to il in the weak topology of measures. Moreover, for all t Cz {to,to + t], we have the 
following convergence: 



limsupsupjlV^^I^™ ■■to<t<to + T} 



hio 
< ess 



sup{|u(t)|„(i) ■.to<t<to + T}. 



(9.1) 



Proof. The only statement which deserves a proof is the last one; if it is not true, 
there exists 7 > 0, a sequence of time steps still denoted by h and a sequence of 
integers m{h) such that 



ymih) 



jj^(h) 



> esssup{|u(i)|Lt) : to < i < ^0 + t} +7- 



(9.2) 



Without loss of generality, we may assume that hm{h) tends to T2 G [0, r]. 

First, T2 cannot be equal to 0: we have learnt in section M that there exists a 
constant C12 and a time ti such that for all h < hi and all m < Ti/h, 

\V"' -V°\ <Ci2mh. 

In particular, this estimate implies that 



ym{h) 



Ijm(h) 



l^oL„ +0{'mh)- 
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but |wol„ is at most equal to {|w(t)|^/j> : ig < ^ < ^o + '''}: which contradicts (9^). In 
the same fashion, we cannot have u{tQ + T2) G int(iir); if it were the case, we could find 
an interval [ti, t^] containing T2 and hi > such that for all h e (0, hi], w/i([ti, t^]) is 
included in a ball of radius r about u(to + T2) included in the interior of K. B ut, i n 
this case, Uh converges uniformly to ii in C*'([ti,T3]) and this contradicts again ( |9.2| ). 
Thus, we assume that T2 is strictly positive and that u{to + T2) belongs to dK. 
Choose a coordinate s ystem such that the origin is at u{to + T2); let ^ be the diffeo- 
morphisni defined at ( 3.13 ). In this case, -D*(0) is given by ( ^.14| ). Define 



Let us compare (3"^ to |l/™|j^,„; it is convenient to define 



F™ = D^{0)C 



then 



We observe that 



|C/™| <|lu,,-u||+C3|m/i-Ti|, 



\X' 



<llt/^ 



and that 



ym _ym 



<C4||dl [2A||X"|| + ||X™-X™-i||]. 



These observations enable us to estimate the difference: there exists a constant C13 
such that 



IT/" 



-/?" 



< Ci3{h + \\u - Uh|lco([fo,io+r]) + I™'* - ^ll)- 



We infer from (3.43) that there exists a constant C14 such that 

/3"+i < min(/3",/3™-i) +Ci4/i. 



We use now (9.2): we can see that for all m < m{h), 

/3"('') < max(/3",/3"-i) + Cii{m{h) - to), 
so that 

max(|y"|f^„ , \V"'^\^_,) > /3™('') - Ci4m{h) - m)h 
Ci3{h + \\u - Uh\\co^ito,to+T]) + I™'* - ^il)- 
If T4 < Ti is such that 

Tl-T4<7/(4Ci4), 
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and if 



Ci3(/l+ \\u - Uh\\co(^[t„^to+r]) + \mh-Ti\) <7/(4Ci4), 



we can see that for all small enough h and all to G { [T4//1] , . . . , m{h)} the following 
estimate holds: 



max(|T/'"|y„,|F'^ 



(9.3) 



> esssup{|u(i)|„(j^ --to <t <tQ + T} +7/2. 
But the function Vh defined by 

Vh{t) = ly"!^™ if t e [mh, (m + l)h) 

converges almost everywhere on [0,t] to \u{t)\u(t)' ^"^ does max(i;/j(f — h),Vh{t)). 
Therefore, in the limit, relation (3.3) leads to 

esssup Vh{t) > esssup{|M(t)|„|-j-| : io < ^ < ^0 + t} +7/2, 

t£lto+T4,to+T2] 

which is a contradiction. D 



A corollary can be inferred imediately from this Proposition and Theorem 3.4: 
Corollary 9.2. For all admissible initial conditions uq and po, there exists 
T > and a solution of 1^, ( pTia] ), (|L4^),(|0]), ( pTSa] ), ( pTSbl ), ( |L8| ) and (|l|) 
defined on [to, to + t]. 

We have proved above the existence of a non-empty interval on which the numer- 
ical scheme converges to a solution of (|1.3[), (1.4a), (1.4b), (1.4c), (1.5a), (1.5b), (1. 



and (1.9). On the other hand, lemma p.l| gives a priori estimates on the solution of 
such a problem. 

We couple now the a priori estimates with the local convergence result to obtain 
a global result: 

Theo rem 9.3. Let R be strictly larger than |po|^ ; O'^d let t{R) be given as 
in lemma S.l. Then, for all small enough h, the solution J7"' of the numerical 
scheme ( |l.l3| ), (1.14), ( 1.16| ), (1.17) is defined on a discrete interval {0, . . . , m(/i)}, 
such that 

hm{h) -^t{R); 

moreover, the approximation Uh converges to a solution u of the continuous time 
equation, i.e. (p^), ( |1.4a| ), ( |l.4b| ),( pr4^ , ( |1.5aD , ( |l.5bD , (|l|) and (pTg]), which is 
defined on [to, to + t{R)]. 

Proof. Let Let {0, . . . , m,{h)} be the discrete time interval for which the numerical 
scheme (1.13), (1.14), (1.16|), (1.17) has a solution; we know from theorem 3.4 that 



Assume that 



liminf /i7Ti(/i) = t > 0. 



T < t{R). 



(9.4) 



It is possible to extract from the sequence {uh)h a subsequence, still denoted by Uh, 
such that on all subinterval [0, r'] included in [O, r] , Uh converges uniformly to u. In 
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particular, thanks to theorem 9.1 we will have 



lim max||y"Mr™ : < m < r/h} 

ft— LI lU 

|2 



and for h small enough we will have 



2 



lim maxJIF""!;,™ : < m < t/K} 

< C2 = 1 + i?max{ M{uy^''^ : \u - uo\ < R} 



Thanks to theorem 3.5, we can find ri such that for all u G K r\B{uo,R), for all t/'~^ 
and U' satisfying condition E{u,ri,C2,h), it is possible to define a solution of the 
scheme for < (m — l)h < r, where r is independent of h. In particular, if we let 

t' = T — t'/2, I = [t' /h\ and u ~ u(t'), we can extend the scheme up to m satisfying 

mh < iiim{lh + t' /2, t{R)), 



which contradicts (9.4). This proves the desired result. D 
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